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/* Example for ratioraetric approach to noise cancelling */ - 

/define LAMBDA 0.95 


/define LAMBDA 0.95 

void OxiLSL_NC( int reset, 

int passes, 

int *signal_l, 

int *signal 2, 

int *signal~3, 

int *target i, 

int *target~2) { 

int if ii* k, m, n, contraction; 
static int *s_a, *s_b, *s_c, *out a, *out c; 
static float Delta_sqr, scale, noise_ref; 

if( reset *= TRUE) { 
s_a ■= signal_l; 
3 s_b = signal_2; 
r= j s_c « signal_3; 
out_a = target_l; 
out_c = target_2; 
r ; factor = 1.5; 
y I scale « l.o I A 160.0; 


* noise canceller initialization at time t=0 */ 

nc[0].berr » 0.0; 
ncfoj .Gamma = l. 0; 

for(m=0; m<NC_CELLS; m++) { 

nc[m] .err_a~ = o.O; 

nc[mj.err_b =0.0; 

nc(m].Roh_a «= 0.0; 

nc[m].Roh_c = o.O; 

nc[mj. Delta = 0.0; 

ncfmj.Fswsqr » 0.00001; 

nc[mj.Bsvsqr «= 0.00001; 


END INITIALIZATION 


for(k=0; k<passes; k++) { 


contraction *= FALSE; 

for(n«0; »< NC_CELLS; m++) { /* update delay elements 

nc[m].berrl « nc[m].berr; 
^ nc[m] .Bswsqrl « nc[m] .Bswsqr ; 


ncise_ref = factor * log(1.0 - (*s_a) * scale) 

- log(1.0 - (*s_b) * scale) ; 

nc[0].err_a « log(1.0 - (*s_b) * scale); 

nc(0].err_b «= log(l.o - (*s_c) * scale); 


# • • • 


++s_a ; 
++s_b ; 
++s_c; 

nc(0] .ferr 
nc(0] . berr 
nc[0] .Fswsqr 
nc[0] .Bswsqr 


noise_ref ; 
noise_ ref ; 

LAMBDA * nc[0]. Fswsqr + noise_ref * noise_rer; 
nc[0] .Fswsqr; 


/* Order Update */ 

for(n»l;( n < NC_CELLS) £& (contraction 

/* Adaptive Lattice Section */ 


FALSE); n++) { 


m « n-1; 
ii- n-1; 

nc[m] .Delta *« 
nc(mj. Delta +« 
Delta_sqr = 


LAMBDA; 

nc[m].berrl * nc[m].ferr / nc[m] .Gamma ; 
nc[m]. Delta * nc(m]. Delta; 


nc[n] . fref 
nc[n j .bref 

nc[n] .f err 
nc(n] .berr 

nc(n] .Fswsqr 
nc[nj .Bswsqr 


= -nc(m]. Delta / nc[m] .Bswsqrl; 
«= -nc[m]. Delta / nc[mj .Fswsqr ; 

= nc(m].ferr + nc[n].fref * nc[m].berrl; 
« nc[m].berrl + nc[n].bref * nc[m].£err; 

■i 

nc[m]'.Fswsqr - Delta_sqr / nc[m] .Bswsqrl ; 
nc(m}. Bswsqrl - Delta_sqr / nc[m] .Fswsqr; 


iff (ncfnl. Fswsqr + nc[n] .Bswsqr) > 0.00001 | | (n < 5) ) { 

nc[n] .Gamma « nc[m] .Gamma - nc[m].berrl * nc[m].berrl / nc[m) .Bswsqrl; 
if (nc[n] .Gamma < 0.05) nc(n]. Gamma « 0.05; 
if (nc[n] .Gamma > 1.00) nc[n]. Gamma =1.00; 

/* Joint Process Estimation Section */ 

nc(m].Roh_a *= LAMBDA; 

nc(m].Roh_a += nc[m].berr * nc[m].err_a / nc[m]. Gamma ; 
nc[m].k_a = nc[m].Roh_a / nc[m) .Bswsqr; 
nc[n].err_a = nc[m].err_a - nc[m].k_a * nc[m].berr; 

nc(m].Roh_c *~ LAMBDA; 

nc[m].Roh_c +«= nc(m].berr * nc(m].err_b / nc[m]. Gamma ; 
nc[m].k_c « nc[m].Roh_c / nc[m) .Bswsqr ; 
nc[n].err_b « nc[m].err_b - nc[m].k_c * nc(m].berr; 

) 

else { 

contraction = TRUE; 
for(i«=n; i<NC_CELLS; i++) { 

0.0; 
0. 0; 
0.0; 
0.0? 
0.0; 

0.00001; 
o.ooooi; 
o . ooooi ; 


nc[i] .err_a 
nc[ i] .Roh_a 
nc[ i] . err_b 
nc[ i] .Roh_c 
ncf i] .Delta 
nc[ i] .Fswsqr 
nc( i) .Bswsqr 
nc ( i ] . Bswsqrl 


) 

} 

} 

*out_a++ = (int) ( (-exp(nc[ii] .err_a) +1.0) / scale) 

*out_c++ = (int) ( (-exp(nc[ii] .err_b) +1.0) / scale) 


) 

) 

/********************** Least Square Lattice 
************************* 


